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Abstract 

We revisit the famous Nose-Hoover system in this paper and show the existence of some averagely 
conservative regions which are filled with an infinite sequence of nested tori. Depending on initial 
conditions, some invariant tori are of trefoil knot type, while the others are of trivial knot type. 
Moreover, we present a variety of interlinked invariant tori whose initial conditions are chosen from 
different averagely conservative regions and give all the interlinking numbers of those interlinked 
tori, showing that this quadratic system possesses so rich dynamic properties. 
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1. Introduction 

In 1984, Shuichi Nose constructed a system called Nose equations to model the interaction of a 
particle with a heat-bath [1,2]. In 1986, Posch, Hoover and Vesely simplified the Nose equations by 
omitting an inessential variable and replaced the residual ” momentum” by a ” friction coefficient” 
and then got the following Nose-Hoover equations [3]: 


J x = y 

< y = —x — yz (1) 

[ z = a(y 2 ~ 1). 

where a is a positive real parameter. System (1) does not possess any fixed points and yet some 
interesting dynamics behaviors have been shown in [3]. For example, for some large value of pa¬ 
rameter a, the system exhibits both tori and chaotic trajectories, depending on initial conditions. 

Recently, Sprott et al. considered generalized Nose-Hoover oscillators (the corresponding har¬ 
monic equations are x = y; y = —x — zy ; z = y 2 — 1 — etanh(x)) [4]. They showed many interesting 
and surprising dynamic behaviors as follows. For some values of e, they got that the conservative 
regions can coexist with dissipative regions in phase space by constructing appropriate cross- 
section. Especially, they showed that the systems owns interlinked invariant tori (see FIG 2 in 
[4]). Also, Sprott have found an other system without equilibrium that has a strange attractor 
and invariant tori [5]. However, it seems that the system studied in [5] does not possesses inter¬ 
linked invariant tori. 

In this paper, we revisit the original Nose-Hoover system (1) with parameter a = 10. By 
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using proper cross-section, we show more types of invariant tori than those shown in [3], such as 
the invariant tori of trefoil knot type, the invariant tori with different types of symmetry, etc. 
Moreover, we show that the original system (1) possesses more complicated interlinked invariant 
tori than the generalized Nose-Hoover oscillators shown in [4] in terms of interlinking numbers 
(see the following Section 3 for definition). 

2. The invariant tori of knot type and periodic orbits in six different averagely con¬ 
servative regions 

Throughout the rest of this paper, we fix a = 10. When taking some initial conditions, the 
corresponding regular trajectories can trace out invariant tori in phase space. For example, for 
this initial condition (x 0 ,y 0 ,z 0 ) = (—0.72,0,0), an invariant torus is produced and its Lyapunov 
exponents are (0,0,0)(see Figure 1(a) and Figure 2). In addition, chaotic trajectories can be found 
for some other initial conditions. For example, when taking the initial condition ( xq,hq,zq ) = 
(—2.4,0, 0), a chaotic trajectory with the Lyapunov exponents (0.0403 -0.0009 -0.0395) is produced 
(see Figure 1(b) and Figure 2). Numerical calculation shows that the chaotic trajectory is a so- 
called attractor. Since the system (1) has no equilibrium points, the attractor is ’’hidden” in 
the sense of Leonov and Kuznetsov [6-8], meaning that its basin does not intersect with small 
neighborhoods of any equilibrium, and it is cannot be found by standard computational methods. 

From Figure 1(c) and Figure 2, it can be seen that the invariant torus and the chaotic attractor 
are coexisting in phase space, which indicates the complexity of invariant sets of system (1). 



(a) An invariant torus 
traced out by initial con¬ 
dition (-0.72,0,0). 



(b) A chaotic attractor 
with initial condition (- 
2.4,0,0). 



(c) The interrelated loca¬ 
tions of the torus in Fig¬ 
ure 1(a) and the attractor 
in Figure 1(b). 


Figure 1: A torus and an attractor of system (1) with a = 10. 

In order to demonstrate more clearly the complicated global structure of invariant sets of 
the system, we construct a cross-section M = {( x,y,z ) e M 3 | y = 0} as shown in Figure 2 and 
study the Poincare map P defined in the plane. Figure 2 shows that the six disjoint subsets 
Di(i = 1,2, ...,6) contained in M are filled with an infinite sequence of nested loops, where D 5 
is the union of two yellow regions D\ and , and D () is the union of two violet regions D\ and 
Hg. Here, each loop is traced out by regular solution with different initial conditions. Moreover, 
we find that, for every initial condition in D t (i = 1,2,3,4), the corresponding trajectory can 
generate a unique loop in Di under the iterated Poincare map . This seems to imply that those 
tori corresponding to the loops from Di(i = 1,2,3,4) are of trivial knot type. The torus shown 
in Figure 1(a) above is just an example of trivial knot type and its initial condition is in D 3 . 
Nevertheless, for each initial condition in Di (i = 5,6), the corresponding trajectory can generate 
two disjoint loops, one of which is in D] and the other is in D~. This should indicate that the tori 
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Figure 2: Detailed cross-section image for a = 10 with 35 initial conditions taken over the interval 
—4.5 < x < 0 with y = z = 0. Here, white regions Tfi,Tb, Dq , Da, , yellow region Dq = D\ U D§, and 
violet region Dq = Dq U Dq are all filled with an infinite sequence of nested loops. 


corresponding to the loops from Dq and Dq are of unusual knot type. For example, the following 
Figure 3 shows two invariant tori of trefoil knot type and their initial conditions are in Dq and 
Dq, respectively. 



(a) An torus of trefoil knot type with 
the initial condition (—1.613,0,0) in 


(b) An torus of trefoil knot type with 
the initial condition (-1,0,-6) in Dq. 


d 5 . 


Figure 3: Two tori of trefoil knot type. 

We note that not all regular trajectories can fill full of the surface of a torus. This is due to a 
well known result in theory of dynamical systems [9], i.e., for the orientation-preserving rotating 
map of circle, if the rotation number is rational, each orbit is periodic and can not £11 full of the 
circle; Otherwise, each orbit can almost fill full of the circle if the rotation number is irrational. 
Thus, for the system (1), considering the Poincare map P restricted to one of the loops in M, 
if the rotation number of this map is a rational number, one regular trajectory produces a finite 
string dots in the loop, which implies that we must choose several appropriate initial conditions 
to produce the closed loop in the cross-section. In contrast, one regular trajectory produces a 
dense set of dots in the loop if the rotation number is irrational. 
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(a) A period orbit with the (b) A period orbit of trefoil 
initial condition in D 4 . knot type with the initial con¬ 

dition in D 5 . 


Figure 4: Two period orbits as examples. 


As we all know, a fixed point of Poincare map corresponds to a periodic orbit of system (1). 
Since the first return Poincare map P is a continuous mapping from D t to Dj(i = 1, 2, 3, 4), there 
exists at least a fixed point in D, (i = 1 , 2, 3, 4) for P by Brouwer fixed point theorem. Neverthe¬ 
less, for i = 5,6, P maps D] to Df and Df to D], Thus, the second return Poincare map P 2 map 
Dj to D\(i = 5, 6), which shows that there exists at least a fixed point in D] (i = 5, 6) for P 2 by 
Brouwer fixed point theorem. From the discussing above, we conclude that the system (1) has 
at least six periodic orbits. In Figure 2, the six small red dots in regions Di, D 2 , D s , Z) 4 , Dl and 
Dl indicate six different periodic orbits of system (1) respectively and two of them are shown in 
Figure 4 above as examples. 

Numerical calculation shows that, if randomly taking xq € D, (i = 1,2, ...,6) , we obtain all 
det(DP(xo)) ~ 1 for x 0 D, L (i — 1,2, 3,4) and det( y DP 2 (xo)) ~ 1 for xq G A (* = 5,6)(See 
Appendix for the Matlab code). This certifies that, for any i = 1,2, 3, 4, 5,6, the phase space 

OO 

domain V) = 1J <pt(D,) which is filled with an infinite sequence of nested tori should obey 

t= — OO 

time-averaged version of the Liouville’s theorem [4], Here, (pt(-) is a flow generated by system 
(1). Thus we treat V) (i — 1,2,..., 6) as averagely conservative regions. 


20 ~i z-axis 



The projection on xy-plane 



(a) An invariant torus contained 
in V 4 which itself is symmetrical 
about z-axis. 





(b) Two invariant tori with inter¬ 
linking number 8 (see Section 3), 
which are symmetrical to each 
other about z-axis. Here the red 
torus is in Vq and the blue torus 
is in V5. 


Figure 5: Two types of symmetry about z-axis. 

Obviously, system (1) is invariant under the transformation of coordinates (x, y , z) -A (—x, — y, z), 
which indicates that, for any invariant torus of system (1), the torus generated by a 180° rotation 
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of this invariant tori about 2 -axis is still an invariant torus of the system. Numerical simulation 
shows that, for system (1), there exists two types of symmetry of tori about z-axis. The first type 
is these tori contained in V)(i = 1,2, 3,4), any of which itself is symmetrical about z-axis, such 
as the torus shown in Figure 5(a) above. The second type is these tori contained in V(5), any 
of which will be transformed into another invariant torus contained in V (6) under the symmetry 
transformation about 2 -axis and vice versa, such as the two tori shown in Figure 5(b). 

In addition, since V • / = —z, all of the tori in V) must go though plane {(x, y, z)\z = 0} (see 
Figure 5). 

3. All interlinking numbers of the different interlinks of tori 

As shown in Subsection 2.1, for any i = 1, 2,.., 6, the tori which stem from D, form an infinite 
nested sequence. However, the interrelated locations of tori derived from different averagely con¬ 
servative regions are still not clear so far. Now, by computer simulation, we show the interesting 
and magnificent interlocked phase-space structures in Figure 6. Here, the six tori shown in Figure 
6 are produced by utilizing the six initial conditions (-2.48,0,0), (-1.97,0,0), (-1.40,0,0), (-0.72,0,0), 
(-1.61,0,0) and (-1,0,-6) which belong to Di, D 2 , D 3 , Z) 4 , D 5 and D 6 respectively, and all their 
Lyapunov exponents are (0,0,0). 





(f) T 23 (g) T 24 (h) T 25 (i) r 26 


(j) T34 



Figure 6: All interlinks of tori derived from six different averagely conservative regions for a = 10. 
Here sub-caption 77, denote that the two interlocked tori in the sub-figure are produced by using initial 
conditions in Di and Dj respectively (1 < i < j < 6). 

From Figure 6, it can be seen that the interlinks of tori exhibit many different types. Following 
[10], we can define interlinking number to measure the complexity for a interlink of two tori 7\ 
and T 2 in R 3 . It can be defined in various ways [11,12], all of which turn out to be equivalent. 
Here, we adopt the easy way used in [10,13] to clarify the idea which is to count up as follows the 
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crossings of Tj and T\ in a ’regular projection’ of the link (a drawing of it such that no more than 
two lines cross at any points). Assign orientations and then put an arrow to the two projections 
of the two tori. Define a number e(p) = ±1 for ’right’ or ’left’-handed crossing p (see Figure 
7(a)).Then define the interlinking number as 

Lk(T u T 2 ) = 

Note that LA:(7j. X 2 ) = Lk(T 2 ,Ti). Figure7(b)[10] shows two examples for the definition of 
interlinking number. 



Figure 7: (a)Definition of the number e(p) for 


crossing p\ (b) Two examples of interlinking numbers. 




(a) Lk{ 1,4) = 2 (b) Lk( 4,5) = 3 


Figure 8: The illustrations for the calculations of LA;(1,4) and LA(1,2). 
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Table 9: All interlinking numbers Lk(i,j) of the interlinks shown in Figure 6 (1 < i ^ j < 6). 

Now, we calculate all the interlinking numbers of the interlinks of tori shown in Figure 6. 
We first stress an obvious fact, namely, for any torus T from D, and any torus T from D 3 (i ^ 
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j ), Lk(T,T) is equivalent. For simplicity, we denote all the same values Lk(T. T) by Lk(i,j). 
According to the above definition of interlinking number, we obtain all thirty values of Lk(i, j)( 1 < 
i 7 ^ j < 6 ) shown in the following Table 9 by choosing applicable projection plane for every 
different interlinks. Two examples are shown in Figure 8 (a) and (b) above for the illustrations of 
the calculations of Lk( 1,4) and LA:( 4 , 5) and both of their projection planes are .xy-planc. 

4. Conclusions 

This article has investigated the Nose-Hoover system and reveals that the co-existence of 
infinite nested tori, a variety of interlinked invariant tori and a chaotic attractor. Furthermore, 
we have provided a mathematical quantitative description for different interlinks of tori from six 
different averagely conservative regions by introducing the concept of interlinking number. It is 
surprising that the quadratic system which seems so simple possesses so abundantly complicated 
dynamic properties. 

It is worth noting that the system ( 1 ), the generalized Nose-Hoover oscillators studied in [4] 
and the unusual system studied in [5] are all three-dimensional quadratic systems and have no 
equilibrium points. Meanwhile, all of them possess invariant tori and chaotic attractors. A natural 
research topic is when the general quadratic systems without equilibrium points possess the co¬ 
existence of invariant tori and chaotic behavior, and moreover possess the interlinked invariant 
tori. In addition, it is obvious that the interlinking numbers shown in Table 8 will be persisted 
under small perturbation of parameter a provided that there are no bifurcations involvement. 
The bifurcation of invariant tori should be another topic to study. 
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Appendix. The Matlab code for the calculation of det(DP(xo)) (The code for det(DP 2 (x o)) 
is similar and we omit it) 
function calculate-D 

options=odeset (’ RelTol ’, 1 e- 9, ’ AbsTol ’, 1 e-10, ’ Events ’, @events); format long 
x0=[-0.72;0;0]; [t,x] =ode45(@rigid 1, [0,le- 6 ] ,x0, []); 

[t, x,te,xe,ie]=ode45(@rigidl,[0,300] ,x(end,:)’, options); 

L=(ie==l); 

PxO=[xe(L,l),xe(L,3)]; 

x01=[x0(l)+0.000001;0;x0(3)]; 

[t,x] = ode45(@rigidl,[0,le-6],x01,[]); 

[t ,x,te,xe,ie] =ode45(@rigidl, [0,300] ,x(end,:) ’,options); 

L=(ie==l); 

Pxl=[xe(L,l),xe(L,3)]; 
x02=[x0(l);0;x0(3)+0.000001]; 

[t,x] = ode45(@rigidl,[0,le-6],x02,[]); 

[t, x,te,xe,ie]=ode45(@rigidl,[0,300] ,x(end,:)’, options); 

L=(ie==l); 

Px2=[xe(L,l),xe(L,3)]; 

Jl(:,l)=Pxl-Px0; 

Jl(:,2)=Px2-PxO; 
det(DPx0)=det( 1000000* Jl); 

function DX = nose(t,x) 

DX =[x(2);-x(l)-x(2)*x(3);(x(2) 2 -l)*10] ; 
function [value,isterminal,direction] = events(t,x) 
value = x(2); 
isterminal = 1; 
direction = 1; 



